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Abstract 

Dark energy accounts for most of the matter-energy content of our universe, yet current theories of its origin rely on radical 
physical assumptions such as the holographic principle or controversial anthropic arguments. We give a better motivated 
explanation for dark energy, claiming that it arises from a small negative scalar-curvature present even in empty spacetime. 
The vacuum has this curvature because spacetime is fundamentally discrete and there are more ways for a discrete 
geometry to have negative curvature than positive. We explicitly compute this effect using a variant of the well known 
dynamical-triangulations (DT) model for quantum gravity. Our model predicts a time-varying non-zero cosmological 
constant with a current value, Ax 10~ 123 in natural units, in agreement with observation. This calculation is made possible 
by a novel characterization of the possible DT action values combined with numerical evidence concerning their 
degeneracies. 
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Introduction 

Multiple independent sets of empirical data [1-4] indicate that 
about 70% of the matter and energy in our universe comes from a 
mysterious repulsive gravitational effect known as "dark energy". 
Understanding the origin of this energy is one of the most 
important problems in physics. Our only current theories involve 
speculative physical assumptions or finely tuned parameters. One 
popular assumption is the holographic principle: the idea that 
the degrees of freedom in a region of space are encoded on the 
region's boundary [5-14]. Other explanations assume the 
existence of exotic matter fields or modify the Lagrangian defining 
general relativity. One recent theory [15] even hypothesizes a 
connection between dark matter and dark energy. See [16,17] for 
reviews of various explanations of dark energy. 

Our work provides a simpler, better motivated model for dark 
energy set within the well-known dynamical triangulations (DT) 
approach to quantum gravity. This model assumes no holographic 
principle, uses no additional matter fields or finely tuned 
parameters, and does not modify general relativity beyond the 
geometric discretization inherent in dynamical-triangulations 
spacetimes. In our model, a positive vacuum energy of the correct 
observed magnitude spontaneously arises from the entropic bias 
toward negative curvature states inherent in DT geometries. Note 
that treating gravity as an emergent mean-field phenomenon 
driven by entropic forces is a popular research perspective at the 
moment [18-24]. 

A reasonable prediction for dark-energy within a quantum- 
gravity theory is only significant if the theory approximates general 
relativity well at large distances. Why should we believe this about 
a theory that uses DT spacetime states? The progenitor of the DT 
theory, called the Regge calculus, has been used successfully in 
numerical general-relativity and quantum gravity for nearly five 
decades [25-31]. The DT model itself [32-36] and its descendent, 
causal dynamical triangulations (CDT) [32,34,35,37-40] have been 



studied for nearly two decades. The numerous successes achieved 
by these theories give confidence that our model can describe 
general relativity at length-scales much larger than Planck's length. 

The model presented in this paper uses the same discretization 
of geometry and the same action as the DT theory. However, it is 
not identical to DT because it puts restrictions on the set of 
triangulations which contribute to the partition function. These 
kind of restrictions are also what distinguish DT from CDT 
although our restrictions are distinct from those in CDT. Note that 
it is not our purpose to advocate "triangulations" as the ultimate 
structure of spacetime. Indeed, in our calculation the discrete 
nature of geometry may be removed at the end without altering 
the predicted vacuum energy. We suspect that the effect described 
in this paper is actually a generic feature of any quantum-gravity 
theory which predicts a discrete spacetime geometry and which 
has general relativity as its large-distance limit. 

Background Material 

General relativity can be written in the Lagrangian formalism 
using the Einstein-Hilbert action, which in natural units is 



1 



[R-2A) +£„ 



Here M is a closed rc-manifold, a Lorentzian metric, R scalar- 
curvature, A the cosmological constant, C m the Lagrangian for 
matter and ^/—gd x the standard volume element. See Table 1 
for a list of commonly used symbols. Note, both R and C m depend 
on g„ v while A does not. Also note that R is the only term in this 
action with a physically distinguished zero value. In quantum field 
theory on a fixed background geometry, an arbitrary constant can 
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Table 1. Meaning of Commonly Used Symbols. 





Symbol 


Meaning 


M 


closed n-manifold 


T 


triangulation of a closed n-manifold 


i 


edge length of all edges in T 


T(M) 


set of all triangulations of M 


Tk(M) 


set of all triangulations of M with K n-simplices 




number of /-simplices in a triangulation T 




average hinge-degree of a triangulation T 




"flat" hinge-degree, fx* n = |r, ^«5.1 (irrational) 


e„ 


dihedral angle in a regular n-simplex, 6„ = cos -1 (l/«) 


A 


cosmological constant 


£«» 


Lorentzian metric 


R 


scalar curvature of g fn , 




Einstein-Hilbert action 




vacuum Einstein-Hilbert action with A = 0 


A R (T,h,... 


Regge action 


A DT (T,l) 


dynamical triangulations (DT) action 


A k 


a DT-action in the TV-action model 


SA 


minimum separation between actions, see equation (14) 


E(T) 


mean DT-action per volume for triangulation T 


E(f>) 


mean DT-action per volume at mean hinge-degree f.i 


E k 


a mean DT-action in the Abaction model 


s t 


spacetime entropy in nats at mean-action 


(£„„„,£•,„,„) 


interval over which mean-actions are regularly spaced 


SE 


minimum separation between mean-actions, see equation (13} 




volume of A'-simplex, all edge-lengths t, Vk(i) = 


V 


spacetime volume V = N n (T)V n (£) 


s* 


Ar-dimensional sphere 


In this table we list some of the commonly used symbols in this paper and their 



meanings. 

doi:1 0.1 371 /journal.pone.0080826.t001 



be added to C„, without changing the observed physics, allowing 
one to simply set A to zero. Thus, it is reasonable to argue, as we 
do in this paper, that the observed non-zero value of A arises from 
quantum effects related to the scalar-curvature field R. 

Hilbert and Einstein showed that the critical points gp, of this 
action satisfy the equations of motion 



(2) 



These are, of course, the field equations for general relativity. 
Here, R^ is the Ricci curvature tensor and the stress-energy 
tensor for matter. In this work we restrict attention to the Einstein- 
Hilbert action for the vacuum with zero cosmological constant 



i 

167t 



R^gcTx. 



(3) 



The critical points of A'^ are metrics which satisfy the vacuum 
field equations. These metrics are Ricci flat everywhere (R^ = 0 at 
every point) and therefore also scalar flat everywhere (R = 0 at 
every point.) Thus these metrics have action exactly zero. Finally, 
in dimensions less than four, the Ricci tensor determines the full 
curvature tensor R m i, so critical points of Ag'^ in these 
dimensions must actually be flat everywhere (R/mX = 0 at each 
point.). 

In his influential 1961 paper [25] Regge proposed a discretized 
version of A E % which applies to triangulated piecewise-linear (PL) 
manifolds. A triangulation T of a closed n-manifold M is a 
combinatorial w-manifold homeomorphic to M given as an 
abstract simplicial complex. Assigning a length £, to each edge 
in T uniquely defines a piecewise-linear metric on T provided 
these lengths satisfy some natural compatibility conditions. If we 
let N/c(T) denote the number of fc-simplices in T, the Regge 
action is given by 



A R (T,£ U ---M(n> 



1 

Urn 



2 ))Vo1(t"- 2 ). 



(4) 



In this equation, the sum runs over all codimension-2 simplices of 
T (called hinges), S(t"~ 2 ) is the total dihedral angle around the 
hinge t"~ 2 , and VoI(t"~ 2 ) is that hinge's volume. It is easy to 
insert a cosmological constant into this action, although here we 
do not. The possibility of incorporating matter fields into Ar is a 
currently active topic of research. See [41—43]. 

Note that Ar has a nice geometric interpretation. The 
summand in this action is the angle deflect in a small triangle 
enclosing and perpendicular to the hinge z"~ 2 eT, weighted by the 
volume of that hinge. Given the close relationship in classical non- 
euclidean geometry between angle defect and curvature, it is 
natural to interpret Ar as a discrete measure of total curvature. 
Because of the success of the Regge action in describing general 
relativity, we will interpret Ar as a discrete measure corresponding 
to the Einstein-Hilbert action, and thus to total scalar- 
curvature. Interpreting Ar as a total curvature is also supported 
by the fact that, like point-wise curvature bounds in Riemannian 
geometry, bounds on the angle-defect for all hinges have profound 
topological consequences for M. See [44-47] for examples. 

The Dynamical-Triangulation Action 

Suppose we fix the abstract simplicial complex T and consider 
Ar as a function of the edge-lengths £ni(T) only. There is a 

large body of numerical evidence [26-29] that the critical points of 
this action define PL-metrics which behave like solutions g flv to the 
vacuum field equations, at least at length scales much larger than 
the maximum edge-length. See [36] for a overview of this work, 
known as the Regge calculus. In this paper, however, we will require 
all edges to have a single fixed length £ so that the action is 
determined only by the structure of T as an abstract simplicial 
complex, i.e. only on the way the simplices in T are attached 
together. This form of the Regge action has been studied 
extensively as part of the dynamical triangulations (DT) approach to 
quantum gravity. We write this action as 



A DT (T/) = A R (T,l, 



. ,£) 



16k 



E (2»- 

n-2 eT 



Meg(T"- 2 )) 



(5) 
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where V^{i) 
of length I, 



l k is the volume of a A>simplex with all edges 



COS _1 (^) is the dihedral angle in such a 
simplex, and clegif), called the degree of x, is the number of n- 
simplices in T with T as a face. Usually, we will suppress the 
dependence on t and write simply Adt(T). 

Now for some terminology and preliminary results. Let T(M) 
denote the set of all triangulations of a fixed closed M-manifold M. 
We will write T k{M) for the set of all triangulations of M 
containing exactly K M-simplices, and Tx,a(M) for those with K 
M-simplices and DT-action A. Since there are only finitely many 
ways to attach together the faces of a finite collection of n- 
simplices, T k(M) and T k,a(M) are finite sets. We define 
S(K,A) = ln|T k,a(M)\ to be the spacetime entropy of M 
for K M-simplices and action A. We will also need notation for the 
average hinge degree of a triangulation T, 



KT) = 



1 



N„- 2 (T) 



E de s(^ 2 )- 



(6) 



By equation (7) we can replace N„-2(T) with ( ~2 I — "t^t an d 



M + 1 

the summation by I 2 )N„(T) to get 



I6n 



Mr) 



w+1 



Finally, moving a factor of 2n ^ ^ )N„( T) to the front finishes 
the derivation. 



Mean Action Per Volume 

The primary observable quantity of concern in this work is the 
mean action per volume, i.e. the average Lagrangian density 
over the manifold: 



By double-counting arguments we may alternately write this as 



N„(T) 
N n - 2 {T) 



N„- 2 {T) ' 



(7) 



Proof. Suppose we examine each (n — 2)-simplex T in T and 
place a mark on each n-simplex with T as a face. Clearly we have 
placed Ylz deg(r) marks. On the other hand, each M-simplex has 
n+l 
2 



codimension-2 faces, so the number of marks is also 



w+1 
2 



N„(T). Dividing through by N n _2{T) gives the first 



equality. Next, suppose we examine each (n — l)-simplex a in T 
and place a mark on each of the two n-simplices incident at a. We 
have obviously placed 2N n -\(T) marks. However, each n-simplex 
has n+l codimension- 1 faces, so the number of marks is also 
(n+l)N„(T) and we have 2N„-\ =(«+ l)N„. Plugging into the 
previous equality and simplifying finishes the proof. 

The first part of equation (7) lets us nicely express Adt(T) as a 
function of the number of n-simplices in T and its average hinge- 
degree. We get 



1 



(8) 



2n 



where ji* n = -— is called the flat hinge-degree. Why do we call 
0„ 

fi* n the flat hinge-degree? It is the number of regular w-simplices 
needed around a hinge to provide a total dihedral angle of exacdy 
271, the expected quantity in a flat space. Note that, except in 
dimension two (where [t 2 = 6) the quantity fi* is not an integer. 

Proof tf Equation (8). We begin with the DT action (5) and 
distribute the sum into the summand to obtain 



A DT (T,£) = 



Vn-im 

1671 



2nN„_ 2 (T)- 



n-2 eT 



(9) 



E(T,£): 



A DT (T,l) 
V(T,i) 



(11) 



where V(Tf)= V n (i)N n (T) is the PL-volume of T. We use the 
symbol E to remind us that this is a physically well-defined global 
observable with dimensions of energy per volume. Equation (8) 
gives a lovely formula for the mean action, 



E(Tf 



1 



1 



(12) 



where C n = 



1 



n+l 
2 



Vn-2(l) 

v„(\) 



depends only on the dimension n. 



This tells us that for a fixed dimension and edge-length the mean- 
action depends only on the average hinge-degree fi. For notational 
convenience we will usually suppress the n and I dependence and 
simply write E(p) or E(T). 

Finally, note that for a fixed number of M-simplices N n we can 
use equations (7) and (12) to find the minimum possible separation 
between mean-actions. This corresponds to changing the number 
of hinges by one, resulting in a change to E of 



8E = D, 



1 



£ n-2 



(13) 



where D„ = - " ^ - and F„ = - K„_2(l) depend only on the 
8 V„(l) 8 

dimension n and V = N„ V„{£) is the total spacetime volume. The 

minimum possible separation between actions is then given by 



5A= V8E = FJ"- 2 . 
Action Spectrum in Dimension Three 



(14) 



From this point forward, we will restrict attention to dimension 
three. What can we say about the possible values of E on Tjf (M) 
when n = 3? This is a formidable problem, since even for a small 
number of tetrahedra K the set T k(M) is quite large and 
complicated. We begin with an elementary result: for any 
triangulation T of a closed 3-manifold M we have 
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N 0 =N 3 --1 



and N\=N 3 - 
H 



(15) 



where Ni = Ni(T) and n = n(T). This means that for a fixed 
number of 3-simplices, the effect of increasing fi (or equivalently, 
decreasing E) is to decrease both the number of vertices iVo and 
the number of edges N\ in the triangulation. 

Proof of Equation (15). We begin with a well-known topological 
fact: every closed 3-manifold has Euler characteristic zero. That is, 
for any triangulation T of a closed 3-manifold M we have 
N 0 (T)-Ni(T) + N 2 (T)-N 3 (T) = 0. Now, we use equation (7) to 
replace N 2 (T) by 2N 3 (T) to get 



N 0 (T)-Ni(T)+N 3 (T) = 0. 



(16) 



so that if N\ >0 is an integer for which E = E(6N 3 /N\) lies in the 
interval (E m i n ,E max ) then E = E{T) for some triangulation T of M 
with JV3 tetrahedra and JVi edges. These E are regularly spaced 
over the entire interval (E m i„,E max ), each separated from the next 
by 



5E-- 



3 1 



1 



(22) 



where V = Ni V 3 (i). This is the smallest possible separation given 
fixed N 3 , so these E are all possible mean-actions on i K E m j n ,E nlax ) . 

Note that in most applications, the number of tetrahedra JV3 will 
be large and the energy densities given in the theorem will be 
approximately 



-o.i9r 



and E m 



iO.ne- 



(23) 



6N 3 (T) 

Using equation (7) again to replace N\(T) by — and then 
rearranging gives 



(17) 



W)=/v 3 (D( — -1 



as desired. Finally, we plug this Nq(T) back into equation (16) and 
simplify to obtain 



Ni(T) = N 3 (T) 



(18) 



completing the proof. 

Equation (15) tells us that to understand the possible values for 
E we must understand the possible combinations of No and N\ 
that can occur in a triangulation of a given closed 3-manifold. A 
1970 paper [48] by Walkup tells us all we need to know. 

Theorem (Walkup). For every closed 3-manifold M there is a smallest 
integer y*(M) so that any two positive integers No and N\ which satisfy 



N) 



>Ni>4N 0 + y*(M) 



(19) 



are given by Ni =N(T) and N 2 = N 2 (T) for some TeT(M). The 
quantity y*(M) is a topological invariant which satisfies 
y*(M)> — 10 for all closed 3-manifolds M. 

Note that y*(M) is known for many manifolds M, see [49], 
although we will not need this information. 

Walkup's Theorem, together with equation (15) and some 
algebra suffice to prove the central mathematical result in this 
paper: 

Theorem. Let M be a closed 3-manifold and N 3 >0 a fixed number of 
tetrahedra. Then, there are mean actions 



N 3 



N 3 -\y\M) 



(20) 



and 



-E\ 6 



N 3 

N 3 +\(i + ^9TW 3 



(21) 



Also note that when edges are Planck's length (£ = 1 in our units) 
the magnitude of these energy densities is enormous, about 10 111 
Joules per cubic meter. 

Proof of Main Theorem. Let M be a closed 3-manifold. W e start by 
showing that if two given integers N\ > 0 and JV 3 > 0 satisfy 



^3+2( 3 + + ^1^3 ( 4 ^3 - f{M)) (24) 

then there is some triangulation T of M with N\ =N\{T) and 
N 3 =N 3 (T). We define N 0 : =Ni -N 3 . Note that N 0 >0 by the 
first inequality in (24). A bit of algebra applied to the second 
inequality in (24) implies 



(25) 



7Y 1 >47Y 0 +7*(M). 



Now, consider the upward opening parabola f{m)- 

m — N 3 which has largest root m 0 = j (3 + % /9 + 8A r 3 ). The first 
inequality in (24) implies Ni — N 3 > \ (3 + V9 + 8iV3) which is just 
No > mo . Since nio is the largest root of an upward opening 
parabola, we conclude /(A f o)>0. By our definition of / and No, 
this tells us 



No 



>Ni. 



(26) 



By Walkup's theorem, inequalities (25) and (26) imply that some 
triangulation TeT(M) has No = No(T) and N\ =N\(T). Finally, 
by equation (15), we know N 3 (T) = N l (T)-No(T) = N 3 as 
desired. 

Next, we divide the inequality (24) by 6N3 and take reciprocals 
to get 



6/V3 



\{AN 3 -f{M)) 



<[i< 



6N 3 



N 3 +U3 + V9 + 8N 3 ) 



(27) 



where n = 6N 3 /N\. Thus, if N 3 is fixed and N\ is an integer for 
which fi = 6N 3 /N\ lies in this interval, then fi = fi(T) for some 
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triangulation T with A^ tetrahedra. By equation ( 1 2) the change in 
mean-action for each increment of N\ is as claimed in equation 
(22), completing the proof. 

The iV-Action Model 

The model used in this paper is designed to be dominated by 
states near a particular chosen target value E* for the mean- 
action. For a fixed number of tetrahedra A^ let Eo be the closest 
attainable mean-action to E* . For each N3, our model admits 
triangulations with mean-action £0 along with those having one of 
the N mean-action values on either side of -Bo- In this paper our 
target mean-action will be E* = 0 since the Einstein-Hilbert action 
for the vacuum in classical general-relativity is zero. Recall that, 
unlike actions in quantum field theory, the values of the Einstein- 
Hilbert and Regge actions are well-defined physical observables. 
This makes such a targeting strategy physically reasonable. 

Why not simply start with a model containing only those 
triangulation T for which E(T) = 0? It turns out that there are no 
such triangulations. That is, for any triangulation T of a closed 3- 
manifold M we have E(T) j^O, or equivalendy ^{T)^ji\. This 
follows from the irrationality of n\ and equation (12). We know ^3 
is irrational due to work [50] by Conway, Radin and Sadun on 
what are called called geodetic angles. Note that these angles are 
actually interesting mathematical objects on their own and are 
central to the solution to Hilbert's third problem on the scissor- 
congruence of polyhedra. 

So, let E-N, ■ ■ ■ ,E// be the mean-actions in the model and 
A-n, ■ ■ ■ ,An the corresponding total actions. Our main theorem 
implies that for any N and spacetime volume V there is an £ small 
enough so that all of the 2N + 1 mean-action values E^ lie within 
the range (E m i n ,E max ) where attainable action-values are regularly 
spaced. For such I our model has partition function 



Z= e s k + i( - A o +k ' 6A > 

k=-N 



(28) 



where Sk = S(N^,Ak) is spacetime entropy at action Ak- The 
expected action for this model is then 



Choosing N 

How are we to choose N? In an ideal world, we would have in 
hand a fully formed DT-style theory of quantum gravity coupled 
to matter, which provably reduced to general-relativity at large 
distances. From this theory we could derive an appropriate N by 
computing how far a typical spacetime was from the classical 
action. We believe such a theory will eventually emerge, but it is 
not yet available. However, we have set up enough machinery to 
reasonably guess what such a theory would tell us about N. 

Suppose we fix a total spacetime volume V and consider the Ab- 
action theory targeting mean-action zero. What happens as we let 
the edge-length £ approach zero? Because the separation between 
actions SA = \£ goes to zero and |„4rj| <&A, if iV is left fixed as 
£— »0 then even the most extreme action values in the theory, 
Ao + N8A, would converge to zero. Since we wish to investigate 
quantum gravity, this is unacceptable and we are forced to choose 
an A^ which diverges as i— >0. Now, suppose we make the affine 
entropy approximation Sk = So + k-rj. Equation (30) implies that if 
rj # 0 then for large enough N and small enough SA the expected 
action is dominated by the final hyperbolic cotangent term and we 
have (A} « sgn(iy)-N<5A This tells us that under these conditions, 
the model is completely dominated by entropy. The oscillating 
complex phase e' -4 which suppresses the contribution of states far 
from A = 0 is swamped by the entropy term involving )j. 

Thus, since 5A is proportional to I, it is natural to choose the 
dimensionless N to be proportional to V 1 ^/£. For such a choice 
we can take the £— >0 limit and the theory gives a finite non-zero 
value for the expected action. Therefore, we choose to use 



N-- 



1 



(31) 



mean-action values on either side of Eq. Notice that by the 
approximations (23) even though A^ diverges as £-*0, all actions in 
the model eventually lie within the "regularly spaced" range 
{E,„i n ,E max } for small enough £. Also note that as £—>0 all the mean 
edge-degrees corresponding to these E, converge to the flat mean 
edge-degree 

Finally, for any fixed rj # 0 we can use equation (30) to compute 
the £->0 limit, obtaining 



<A}=- J2 {A 0 +k-SA)e s k +i(A o +kdA \ 



(29) 



k=-N 



A Euclidean version (A} Euc °f this expected value can be found 
by applying the standard Wick rotation i \— > — 1 to the expression 
above. 

It is currently impossible to write (A} or <_4> Euc as exact 
closed form expressions since the entropies Sk are beyond our 
ability to compute. However, if we replace Sk with its first order 
approximation Sk = So+k-ri for r\ a constant, then a closed-form 
expression can be found. We used the computer-algebra package 
Mathematica to show 



<A} = A 0 



6A SA 
e n + i&A__\ + e {2N+\)(n + iSA) _ \ 

+NdAco\h[(2N+ Y)(n + idA)\. 



(30) 



A closed form expression for <»4>euc can obtained as before by 
replacing i with — 1 in the equation above. 



(32) 



For r] = 0 we get a purely imaginary standard expectation 



lim^o (Ay = i 



1 1 \ 1 

l--F*cot -K5 



and a Euclidean expecta- 



tion given by lim^_>o <-4> euc = 1 — o ^' cot ^ ' " V * 



Evidence for the Entropies St 

The calculation of the expected action as £— >0 given by 
equation (32) depends on two assumptions about the entropies Sk- 
First, for the states contributing to the model, spacetime entropy 
must be an approximately linear function of mean-action, i.e. 
S/cXSo + k-rj, at least for large enough A^. Second, this rj must not 
approach zero as A^— >go. In this section, we present evidence 
from Monte-Carlo simulations and small-A^ enumerations that 
strongly supports these assumptions. 
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Monte-Carlo Sampling Results 

To measure the dependence of entropy on mean-action we use 
a Metropolis-Hastings algorithm to take samples T from T(M) 
near a given number of tetrahedra and mean-action. The 
algorithm wanders among the elements of T(M) by using the 
well-known Pachner moves to change from one triangulation to 
another, repeatedly choosing a random move and executing it 
with probability e~ AU where U is some non-negative objective 
function. Metropolis proved that if we wait long enough between 
samples, then each sample T occurs with probability e~ u ^ T \ 
Here, we use a quadratic objective function 



i2u(E-E 1: )5E. 



(34) 



U( T) = a(E( T)-E r ) 2 + P (n 3 ( T) - 



■Nl 



(33) 



with a>0 and /?>0 fixed constants. This form for U keeps the 
sampled triangulations near a target mean-action and number 
of tetrahedra iV-j . 

If there were equally many triangulations at each E and N3 then 
our sampled pairs (E,Ni) would form a Gaussian distribution 
centered at the target point (E^ ,Nt). If our samples have a 
Gaussian distribution but with mean (E,N3) significantly displaced 
from the target, this indicates a linear dependence of spacetime 
entropy on E and N-j with the magnitude of the dependence 
proportional to the size of this displacement. Since it is obvious 
that spacetime entropy is strongly dependent on Nj and because 
the relative deviation from the mean for N3 is at most « 1 % in our 
data, we focus solely on deviation in mean-action E — E^. From 
this we can estimate the change entropy (per mean-action step) r\, 
in nats, using 

700 



The sampling trials conducted for this paper use the 3-sphere 
M = 5 3 with target mean-action zero (E^ = 0) and various targets 
for the number of tetrahedra 500 < < 7000. In all cases, we 
take a = 3.5xl0 6 and /?=1.0xl0~ 2 . In order to ensure 
independent samples, the algorithm attempts Pachner moves until 
x 10 accepted moves per tetrahedron have occurred. We checked 
that this wait time was sufficient using standard correlation tests. 
For these parameters, each sample was uncorrelated from the 
next. We also checked that the sampled JV3 and E were 
independent. As desired, samples are approximately normally 
distributed with sample mean E somewhat displaced from the 
target £' =0. This indicates that entropy is approximately a linear 
function of mean-action near E = 0, as was assumed in the 
previous section. See Figure 1 for a histogram of mean-actions for 
2700 samples at TV] = 1 70 1 . For each such distribution, we use 
equation (34) to infer the approximate change in entropy r\ 
between mean-actions. These r\ are comfortably negative and do 
not appear to approach zero as N\ gets larger, validating our 
second assumption. See Figure 2. Copies of the code used for 
triangulation sampling are available on request. 

Triangulation Census Data 

In addition to Monte-Carlo sampling evidence, one can also 
see a bias towards negative action states in computer-generated 
censuses of 3-manifolds triangulations. In particular, recent 
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Figure 1. Monte-Carlo sampling of triangulations of V near mean-action zero. We plot the distribution of mean actions E at 1= 1 for 2700 
sampled triangulations of the 3-sphere S 3 . Samples were obtained from a Metropolis-Hastings algorithm using Pachner moves and a quadratic 
objective function «(£-£ t ) 2 + /i(A' 3 -A'|) 2 targeting £ + = 0 and A?f = 1701 with a = 3.5 x 10 6 and fj= 1.0 x 10~ 2 . Waiting times were chosen so that 
Rj 10 accepted moves per tetrahedra occurred between successive samples. Observed means were ^3 = 1770 with standard deviation <7 = 7.0 and 
E= —6.3 x 10~ 4 with standard deviation <r = 3.7 x 10~ 4 . Note that I and E are given in Planck units, L p and E p j L? p respectively. 
doi:1 0.1 371 /journal.pone.0080826.g001 
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Figure 2. Entropy remains a decreasing function of mean-action as the number of tetrahedra grows. We plot the change in spacetime 
entropy, ;/ in bits, due to each minimal increase SE in mean-action for the 3-sphere S 3 near £ = 0, versus mean number of tetrahedra JV3. Values were 
inferred from the bias seen in Monte-Carlo samples of triangulations near £ = 0. See Figure 1 . All data points except the last two were computed from 
2700 samples. At the two largest 7V 3 values, we used 2394 and 1108 samples respectively. Error bars indicate 95% confidence intervals. 
doi:1 0.1 371 /journal.pone.0080826.g002 



advances in enumeration algorithms have allowed for the 
creation of an explicit list of all triangulations of any closed 3- 
manifold using at most 11 tetrahedra. See [51,52]. Unfortu- 
nately, the definition of a "triangulation" used in these censuses 
is slightly more general than ours. They define a triangulation 



of a closed 3-manifold M as a space homeomorphic to M 
obtained by identifying the faces of some finite set of 
tetrahedra. We believe this to be a largely technical distinction, 
and we expect this data to provide a good guide to the general 
features of our set of triangulations T(M). See Figure 3 for a 




°-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Mean Action E in E p /L* 

Figure 3. Entropy versus mean-action from triangulation census data. We plot spacetime entropy S in bits for the three-sphere S* 3 at various 
numbers of tetrahedra YV3, versus mean action £ at 1= 1. Data come from a complete census [51,52] of the s;47 million triangulations of S 3 with at 
most 9 tetrahedra. Note that I and E are given in Planck units, L p and E p /L p respectively. 
doi:1 0.1 371 /journal.pone.0080826.g003 
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graph of spacetime entropy versus mean-action for the 3-sphere 
S 3 and 5<Ni<9. We observe two trends in the data. First, as 
we expect, the number of triangulations increases as the 
number of 3-simplices grows. However we also see the same 
effect as observed in the Monte-Carlo sampling experiments: 
the number of triangulations at a given action is a decreasing 
function of action. 

The Origin of Dark Energy 

Taking rj < 0 and dividing through by V in equation (32) gives 
Jim<£>=!im<£> Euc = -iH. (35) 



Thinking of them as having explanatory power on their own is 
surely deserving of the label "numerology". However, this epithet 
should not be applied to a physically well-motivated theory which 
predicts ab-initio such a numerical relationship, as our model does. 

Discussion 

Our derivation of A has some interesting features. Using the 
Hubble parameter to define our characteristic volume V means 
that the model actually predicts a time-varying cosmological 
constant 

A(t) K \-H(ti (38) 



Let us briefly discuss the physical meaning of Our goal 

was to construct a theory dominated by states close to the classical 
value of the mean-action, E = 0. We did this by "slicing" the 
partition function according to action-value, retaining only states 
whose actions lie within a certain distance of zero. If the volume of 
spacetime is large compared to Planck's volume then we come 
very close to accomplishing our goal. That is, for V » 1 we do 
indeed obtain (E~)xO in the I— >0 limit. However, there is a small 
perturbation away from zero because of the relative entropy of 
action values. Notice that since action values are global 
observables, this effect is independent of the local details of the 
"metric", i.e. the local structure of the triangulation. This leads us 
to expect that, for a typical triangulation at a given £, the average 
action will appear very uniform at length-scales much larger than 
I. Finally, recall that everything in the Einstein-Hilbert action 
except the cosmological constant A depends on the metric g^. 
Thus, the basic structure of Aeh almost demands we interpret our 
non-zero <£) as an emergent cosmological constant given by 

A=-\(Ey=^V-i (36) 

We now turn to the question of applying this result to our own 
universe. This is a somewhat speculative endeavor since our world 
appears to be both 4-dimensional and infinite in extent. However, 
as an entropic effect connected with the pattern of attachment 
between simplices, we expect the perturbation away from E = 0 
identified in this paper to occur quite generally. So, what V is 
appropriate for assessing the magnitude of this effect in our 
particular universe? Considerations of causality give us a 
reasonable answer: take the volume of space which has had time 
to causally communicate with our point of observation. That is, we 
ought to use something like the current Hubble volume H^ 3 where 
Hq is the Hubble constant. Plugging in H 0 x\.2x 10" 61 in Planck 
units gives 

/i«i(r 123 ( 37) 

which is in general agreement with observation. 

At this point, we feel obliged to briefly discuss the term 
"numerology". It has long been known that the observed 
cosmological constant was approximately Hq. This and many 
other unexplained approximate numerical relationships between 
cosmological parameters are often called large number coincidences. 



where H(t) is the Hubble parameter at proper time t. That is, 
we predict that A scales like the area of the cosmic horizon. 
Amazingly, although we made no holographic assumption, this is 
the same behavior that emerges from holographic dark energy (HDE) 
theories [5,6,10,12,13]. In fact, our model shares several other key 
features with these approaches, including the presence of two "cut- 
offs" in the theory which are removed in a coordinated fashion. 
HDE models typically contain both a UV and IR field cut-off 
which are removed in a way that saturates entropy in the 
Bekenstein bound. In our theory, the cut-offs I and N are chosen 
to keep the entropic perturbation on <£> bounded as £-*0. While 
HDE theories are very different in detail from our model, the 
broad similarities are quite striking. Perhaps both approaches are 
pointing to the same underlying physical issues. We hope that the 
relative simplicity of our model can help elucidate these issues. 

We should also mention another explanation for A which shares 
some features with our approach. In [53] it is argued that the true 
ground-state vacuum has A = 0 but that we observe A > 0 because 
the universe has not yet had time to decay into this ground state. 
The author considers a model in which the true ground state is 
given by the superposition of two degenerate A > 0 states, one of 
which describes the universe's present-day vacuum. Since the 
decay probability in a given volume and time period is related to 
the energy density A, the requirement that no such decay has yet 
happened in the Hubble volume provides an estimate for A which 
agrees with observation. This argument leads, as does our model, 
to a connection between the Hubble parameter and A. Also note 
that both models contain states at or near A = 0 which are 
suppressed compared to the A > 0 states. 

Finally, we note that in the very early universe our model 
predicts large A and hence rapid expansion. This raises the 
tantalizing possibility that big-bang inflation and dark-energy are 
manifestations of a common effect, though it is likely that a more 
sophisticated choice for the characteristic volume V would be 
needed. See [54] for consideration of this idea in the HDE context. 
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